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The  following  section  sumtoarizes  the  projects  supported  by  the  contract. 

A.  An  approximate  method  for  calculating  the  reflection  and  transmission  coef- 
ficients of  obliquely  incident  Rayleigh  waves  across  two  welded  quarter-spaces 
is  outlined.  The  transmitted  energy  follows  a reciprocity  relation  and  de- 
creases rapidly  away  from  normal  incidence;  the  reflected  energy  is  small  and 
decreases  with  increasing  angle  of  incidence.  This  numerical  method  is  very 
fast  and  has  a low  memory-storage  requirement. 

B.  A moving-window  analysis  was  employed  to  study  the  group  velocities  of  funda- 
mental mode  Rayleigh  and  Love  waves  with  paths  traversing  the  Tibetan  Plateau. 
The  models  derived  from  best-fitting  these  observed  dispersion  curves  indicate 
that  the  Tibetan  Plateau  has  an  unusual  crust  about  70  km  in  thickness,  with 
a low  velocity  zone  at  an  intermediate  depth  within  the  crust. 


C.  Travel-time  residuals  determined  from  records  of  the  L-DGO  seismic  network  in 
the  northeastern  U.S.  were  used  to  map  the  extent  of  a suspected  suture  zone 
in  northern  New  York,  central  Vermont,  and  southwestern  Quebec.  The  results 
indicate  that  a zone  of  large  negative  residuals  exists  in  northern  New  York 
State.  Recent  data  also  show  an  azimuthal  dependence  in  the  residual  pattern. 

D.  Fault-plane  solutions  of  moderate-sized  earthquakes  along  the  North  Atlantic 
Ridge  were  determined  from  surface-wave  amplitude  spectra  constrained  by  body- 
wave  first  motions.  The  strikes  of  the  faults  derived  from  surface  wave  studies 
agree  well  with  the  local  orientation  of  the  median  valley  near  the  epicenters. 
The  study  demonstrated  that  amplitudes  can  be  used  successfully  to  obtain  re- 
liable focal  mechanism  solutions  for  oceanic  paths  for  periods  from  25  to  80  sec. 

I 


reveals  that  the  French  test  site  In  Algebra  has,  contrary  to  previous  belief, 
not  been  tectonically  stable  since  the  Frecambrlan.  The  site  Is  located  In 
a zone  that  has  experienced  tectonic  activity.  Including  volcanlsm  and  uplift, 
during  the  past  25  million  years.  Seismic  waves  leaving  this  region  may  be 
more  attenuated  In  amplitude;  consequently,  yield  estimations  for  the  French 
tests  will  have  to  take  this  new  finding  Into  account.  The  attenuation  of 
waves  from  that  test  site  may  be  more  similar  to  that  of  waves  leaving  the 
Nevada  Test  Site. 

The  following  section  summarizes  projects,  not  directly  supported  by  the  contract, 

but  of  Interest  to  It: 

A.  The  complicated  tectonics  of  central  Asia  have  been  studied  In  several  aspects. 
The  subjects  covered  under  this  research  Include: 

1.  Detailed  modelling  of  the  Hazara  region  of  northern  Pakistan  based  on  the 
mapped  distribution  of  seismicity,  composite  fault-plane  solutions,  and 
geological  evidences. 

2.  Large-scale  selsmotectonlc  trends  for  a portion  of  south-central  Asia  based 
on  a compilation  and  critical  review  of  the  historical  and  modem  seismicity. 

3.  Interpretation  of  the  Makran  region  of  southern  Pakistan  as  a zone  of  active 
subductlon  from  a detailed  comparison  of  tectonic  features  In  the  Makran 
with  those  In  a subductlon  environment. 

4.  A temporal  and  spatial  study  of  seismic  activity  prior  to  the  1945  earth- 
quake (Ms  ■ 8)  along  the  Makran  coast  of  southern  Pakistan. 

B.  Several  studies  concerning  the  seismicity  and  tectonics  of  the  Aleutian  Islands 
were  carried  out  at  L-DGO.  The  results  of  these  studies  may  be  useful  In  under- 
standing other  subductlon  zones  such  as  the  one  along  the  Kurile  Islands  and 
Kamchatka. 

1.  Hypocenters,  determined  from  data  recorded  by  a local  network  In  the  Shumagln 
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Islands  Alaska,  define  a Benloff  zone  that  Is  about  10  km  thick  and  extends  ^ 

to  about  180  km  depth.  Large  earthquakes  appear  to  rupture  that  shallow-  I 

dipping  portion  of  the  subductlon  zone  that  occurs  above  40  km  depth;  below  | 

that  depth  the  Benloff  zone  is  congruent  along  the  entire  Aleutian  arc.  This 
observation,  and  the  low  level  of  seismic  activity  at  the  present  time  in  the 
main  thrust  zone  compared  to  that  in  the  Benloff  zone  below  40  km,  suggest 
that  seismic  activity  in  the  main  thrust  zone  is  dominated  by  great  earth- 
quakes and  their  aftershocks,  followed  by  long  intervals  of  quiescense. 

2.  Two  high  stress-drop  earthquakes  from  the  Shumagln  Island,  Alaska  were  studied 
in  detail  using  data  from  the  local  short-period  stations,  strong-motion 
accelerograph  (SMA),  and  the  long-period  WWSSN  stations.  The  distribution 

of  the  malnshocks  and  associated  aftershocks  were  located  and  mapped  with  j 

high  precision.  The  SMA  data  were  anlyzed  to  give  detailed  source  models.  ' 

A tectonic  interpretation  for  these  earthquakes  is  given.  The  inferred 

I 

stress  drops,  about  500  bars,  are  among  the  highest  values  known  to  be 
associated  with  any  earthquake.  ! 

3.  Seismicity  in  the  Adak  region,  Alaska  shows  two  thin  zones  about  25  km  apart 

below  a depth  of  100  km.  These  zones  merge  at  a depth  of  175  km.  Observed  i 

focal  mechanisms  showing  down-dip  compression  in  the  upper  zone  and  down-  | 

dip  tension  in  the  lower  zone  are  predicted  as  the  result  of  unbending  of  * 

an  elastic-plastic  lithosphere.  I 


C.  A focal  mechanism  study  and  aftershocks  survey  were  undertaken  to  investigate 
a large  intraplate  earthquake  about  350  km  southwest  of  Bentuda.  Preliminary 
results  are  reported.  That  shock  was  situated  on  old  oceanic  crust  where  the 
attenuation  of  seismic  waves  is  likely  to  be  small.  Large  high  frequency  Fn 
and  Sn  waves  were  recorded  in  eastern  North  America  at  distances  of  about  15* 


to  20 
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D.  A comparative  surface-wave  technique  was  used  to  precisely  determine  the  depths 
of  two  earthquakes  that  occurred  seaward  of  the  Tonga  trench.  A focal  depth  of 
49  km  was  found  for  a thrust- faulting  earthquake,  and  a depth  of  14  km  for  a 
normal  faulting  earthquake.  The  Implied  stresses  are  consistent  with  bending 
of  an  elastic  plate  model. 

1 
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INTRODUCTION 


► 
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During  the  past  two  years  the  Advanced  Research  Projects  Agency  (ARPA)  has 
supported  research  at  the  Lamont-Doherty  Geological  Observatory  of  Columbia  Uni- 
versity under  a contract  entitled  "Suture  zones,  seismic  wave  propagation,  and 
the  tectonics  of  central  Asia".  Research  was  undertaken  to  Investigate  the  fol- 
lowing specific  problems: 

1.  The  crustal  and  upper  mantle  structure  beneath  the  Tibetan  plateau, 

2.  The  effect  of  a vertical  discontinuity  on  the  propagation  of  obliquely  Inci- 
dent Rayleigh  waves, 

3.  The  characteristics  and  size  of  zones  with  anomalous  travel-time  residuals 
that  may  be  caused  by  ancient  suture  zones, 

4.  The  focal  mechanisms  of  small-magnitude  earthquakes  (m^  < 5.S)  determined  by 
an  analysis  of  surface  waves, 

5.  The  tectonic  setting  of  the  Algerian  test  site. 

In  addition,  several  other  studies  have  been  carried  out  at  Lamont  that,  although 
supported  by  other  contracts,  are  of  direct  Interest  to  the  subject  of  this  con- 
tract. These  studies  Include  detailed  work  on  the  tectonics  of  a specific  area 
In  central  Asia,  regional  studies  of  central  Aslan  tectonics.  Investigations  of  seis- 
micity In  the  Aleutians,  preliminary  results  concerning  an  Intraplate  earthquake 
southwest  of  Bermuda,  and  the  precise  determination  of  depth,  employing  surface-wave 
techniques,  for  two  earthquakes  seaward  of  the  Tonga  trench. 

The  research  program  outlined  above  has  Improved  our  knowledge  on  the  dla- 

1 

crimination  of  nuclear  explosions  from  natural  earthquakes  In  several  aspects.  j 

For  the  source  region,  the  program  has  clarlflad  our  understanding  on  the  tectonic  1 

envlronsMnt  In  which  several  nuclear  explosions  were  detonated.  For  the  propagation  | 

path.  It  has  attempted  to  quantify  the  effects  various  tectonic  structures  may  have 
on  the  travel-time  of  body  waves  and  aBq>lltude  of  surface  waves.  The  program  has 

4 
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also  extended  the  usefullness  of  an  established  discrimination  method  to  events  I 

with  small  signal  to  noise  ratio.  In  the  following  sections  the  methodology  and  ^ 

results  of  the  individual  research  projects  are  discussed.  ! 


♦ 


A.  Reflection  and  Transmission  of  Obliquely  Incident  Rayleigh  Waves  at  a Vertical 
Discontinuity  Between  Two  Welded  Quarter-Spaces 


An  approximate  method  for  calculating  the  reflection  and  transmission  coef- 
ficients of  normally  incident  Rayleigh  waves  across  two  welded  quarter-spaces 
was  previously  outlined  by  Alsop,  Goodman,  and  Gregersen  (1974).  Their  method 
exploits  the  orthogonality  property  of  eigenfunctions  in  a layered  half-space 
to  calculate  the  degree  of  matching  between  an  eigenfunction  and  a wave  field 
(Herrera,  1964).  Chen  and  Alsop  (1979)  have  revised  two  of  the  concepts  presented 
in  Alsop  et  al.  (1974),  then  extended  the  revised  method  to  allow  for  oblique 
incidence.  Results  for  normal  incidence  are  in  agreement  with  previous  experi- 
mental and  theoretical  results  (McGarr  and  Alsop,  1969;  Alsop  et  al. . 1974).  This 
method  features  two  definite  advantages,  speed  and  low-memory  storage  requirement, 
over  other  numerical  methods  commonly  used  on  surface  wave  transmission  problems 
such  as  finite-element  or  finite-difference  methods. 

For  two  media  with  a phase  velocity  ratio  of  1.16,  they  found  that  the  trans- 
mitted energy  follows  a reciprocity  relation;  that  is,  if  the  same  propagation  path 
is  followed,  the  energy  transmitted  from  medium  I to  medium  II  is  the  same  as  that 
from  medium  II  to  medium  I.  Theoretically,  this  reciprocity  relation  should  exist 
between  any  discrete  eigenfunction  on  either  side  of  a transition  zone.  The  trans- 
mitted energy  decreases  from  near  lOOZ  at  normal  incidence  to  50%  at  about  40” 

(Figure  la).  The  reflected  energy  is  less  than  1%  and  decreases  with  increasing 
angle  of  incidence  (Figure  lb).  Extending  their  calculations  to  refracted  angles 
beyond  the  critical  angles  of  S to  P and  S to  S conversions,  they  found  that  the 
transmitted  energy  continues  to  follow  the  reciprocity  relation  and  decreases  further. 
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Reflected  energy  remains  below  IZ  but  fluctuates  Irregularly;  no  total  internal 
reflection  Is  observed. 


Since  Chen  and  Alsop  did  not  take  boundary  conditions  on  the  free  surface  Into 
account,  they  believe  that  the  calculated  energy  for  reflection/transmission  will 
be  reduced  by  diffractive  effects,  and  their  results  probably  represent  an  over- 
estimate for  the  reflected  and  transmitted  energy.  They  are  presently  working  on 
the  calculation  of  Rayleigh  waves  across  a vertical  discontinuity  In  two  layered 
quarter-spaces,  and  will  compare  the  theoretical  results  with  observations  from 
an  ocean-bottom  seismograph  and  an  adjacent  land  station. 

B.  Crustal  and  Upper  Mantle  Structure  of  the  Tibetan  Plateau 

A knowledge  of  the  velocity  and  attenuation  structure  beneath  the  Tibetan 
Plateau  Is  Important  for  the  discrimination  of  explosions  detonated  In  the  Sln- 
klang  Province,  China  and  eastern  Kazakh,  U.S.S.R.  To  determine  crustal  models 
for  the  Tibetan  plateau  Chun  and  Yoshll  (1977)  and  Chun  (personal  communication, 

1977)  determined  the  group  velocities  of  fundamental  mode  Rayleigh  and  Love  waves 
In  the  period  range  5-100  sec.  A moving  window  analysis  was  employed. 

The  models  derived  from  these  observed  dispersion  curves  (Figure  2)  Indicate 
that  the  Tibetan  plateau  has  an  unusual  crust  about  70  km  In  thickness.  Crustal 
velocities  are  generally  low  (Figure  3 ) . From  the  abundant  volcanoes  in  the  re- 
gion, the  low  velocities  are  Inferred  to  be  caused  by  high  temperatures  at  crustal 

I 

depth.  A low  velocity  zone  at  an  Intermediate  depth  within  the  crust  appears  to  | 

be  well  constrained  by  the  observed  dispersion  data. 

C.  Travel-T<n»*  Study  of  Anomalous  Crust  and  Upper  Mantle  Structure  In  the  North- 
eastern United  States 
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Slabs  of  mantle-llke  material  Imbricated  In  the  upper  crust  during  the  suturing 


t 

I 

of  two  continental  masses  can  result  in  an  anomalous  zone  of  high  velocity  mat-  | 

erial.  Seismic  waves  that  pass  through  such  zones  will  be  affected  in  travel- 
time, and  possibly  in  waveform  and  amplitude.  The  delineation  of  their  effects  j 

are  important  for  Che  accurate  locating  of  nuclear  explosions  and  the  computation 
of  their  yields. 

Travel-time  residuals  determined  from  records  of  Che  L-DGO  seismic  network  in 
the  northeastern  U.S.  were  used  to  map  the  extent  of  a suspected  suture  zone  in 

northern  New  York,  central  Vermont  and  southwestern  Quebec.  The  travel-time  real-  j 

j 

duals  are  computed  for  Impulsive  F-wave  arrivals  from  both  teleselsmlc  earthquakes 
and  explosions.  Preliminary  residuals  are  computed  for  each  station  by  using  the 
Jeffreys-Bullen  tables.  The  average  of  these  preliminary  station  residuals  is 

,1 

Chen  subtracted  from  each  individual  residual  to  produce  the  final  residual.  In 

this  way  errors  caused  by  source,  path  or  other  effects  common  to  all  stations,  ! 

hopefully  are  removed.  , 

To  date  eighteen  events  have  been  examined.  The  results  tend  to  confirm  the 
observations  of  Fletcher  ^ al.  (1978)  Chat  a zone  of  large  negative  residuals 
exists  in  northern  New  York  State  and  along  the  St.  Lawrence  River  (Figure  4).  The 
recent  data  indicate  that  the  residual  pattern  is  azlmuthally  dependent.  The  reel-  ’ . 

duals  range  from  -1.8  to  -0.5  sec  for  events  to  the  north  and  west,  but  are  negll- 

I 

glble  for  other  directions.  Data  previously  suggestive  of  a similar  zone  of  large  t 

r 

negative  residuals  in  the  greater  New  York  City  region  are  conflicting  in  the  light  b 

of  recent  data.  i 

( 

The  zone  of  negative  residuals  in  northern  New  York  State,  central  Vermont,  and  | 

I I'. 

southern  Quebec  may  result  from  mantle  material  chat  was  thrust  into  the  crust  in  j 'f ' 

this  region  during  the  suturing  of  the  rest  of  North  America  and  New  England  in  Or-  | ; 

dovlclan  times.  Outcrops  of  mafic  material  in  this  region  support  such  an  inter-  | 

i I- 

pretatlon.  Further  research  is  still  needed  to  determine  if  anomalous  amplitudes  N 

f ! 

and  waveforms  can  be  documented  and  to  apply  this  method  to  other  regions  where  su- 

Cure  zones  are  suspected  (such  as  the  numerous  suture  zones  of  Asia) . 1 


D,  Surface  Wave  Studies  of  Small-Magnitude  Earthquakes  from  the  Mid-Atlantic  Ridge 

For  earthquakes  from  the  mid-ocean  ridges,  a combination  of  conditions  at  the 
source  and  receiver  often  makes  focal  mechanism  studies  using  first  motion  dif- 
ficult. At  the  source,  the  moderate  size  of  the  mid-ocean  ridge  earthquakes  and 
the  strong  attenuation  of  the  low-velocity  zone  there  combine  to  result  In  small 
and  emergent  first  motions  (Solomon,  l'^73).  In  addition,  most  long-period  obser- 
vational stations  are  located  at  teleselsmlc  distances  from  mid-ocean  ridges,  so 
that  most  stations  fall  Into  a dllatlonal  quadrant  for  normal-faulting  events.  Thus, 
the  orientations  of  the  nodal  planes  are  often  difficult  to  resolve. 

Recent  developments  In  normal  mode  theory  (Salto,  1967)  and  signal  processing 
(e.g.  Landlsman  ^ al . , 1969)  have  made  the  analysis  of  surface  waves,  which  usually 
have  the  largest  amplitude  on  a long-period  seismogram.  Into  a more  common  method  In 
the  study  of  fault-plane  solutions  (Forsyth,  1975).  There  Is,  however,  a major  draw- 
back In  the  surface  wave  approach  In  that  it  requires  a large  amoun  of  signal-proces- 
sing. Mid-ocean  earthquakes,  by  nature  of  their  peculiar  tectonic  setting,  require  less 
data-processlng  when  using  the  surface  wave  approach  to  obtain  fault-plane  solutions. 

It  Is  well  knoxm  that  mid-ocean  ridge  earthquakes  have  either  normal-  or  strlke- 
sllp-faultlng  mechanisms  (Sykes,  1967),  and  they  are  usually  shallow  (Weldner  and 
Aki,  1975).  Normal  mode  theory  for  Rayleigh  waves  predicts  that  for  shallow  events, 
normal-faulting  will  generate  a two-lobed  radiation  pattern  with  amplitude  nodes 
along  the  fault  strike,  but  strlke-sllp-f suiting  will  generate  a four-lobed  radia- 
tion pattern  with  amplitude  nodes  along  the  strikes  of  both  conjugate  nodal  planes. 
Therefore,  Rayleigh  wave  radiation  patterns  from  mid-ocean  earthquakes  will  be 
Informative  to  their  fault-type  and  approximate  orientation  of  fault  strike.  Ho 
and  Sykes  have  used  this  property  to  study  the  focal  mechanisms  of  some  moderate- 
sized (4.6  Ms  £ 5.6)  earthquakes  In  the  North  Atlantic  Ocean  (Table  1,  Figure  5A) . 


Long-period,  vertical  component  seismograms  from  the  WWSSN  were  selected  from 
circum-Atlantic  stations.  The  seismograms  were  then  digitized  and  interpolated  at 
1-s  intervals.  A moving-window  analysis  (Landismen  ^ • , 1969)  with  a window 

length  7 times  the  period  was  performed  on  each  record  to  obtain  the  amplitude  spec- 
tra. They  then  corrected  the  amplitude  spectra  for  Instrument  response,  geometrical 
spreading  and  attenuation  using  a general  Q model  (Tsai  and  Akl,  1969).  The  spectral 
densities  at  any  specific  period  were  then  plotted  as  a function  of  eplcentral 
azimuths  (Figure  5B-E) . The  amplitude  nodes  in  the  Rayleigh  wave  radiation  pattern, 
together  with  constraints  from  body  wave  first  motions,  allow  them  to  resolve  the 
fault-type  and  the  strike  of  the  faults  to  + 15®.  The  strikes  of  the  faults  derived 
from  surface  wave  studies  agree  well  with  the  orientation  of  either  the  medlam 
valley  or  fracture  zones  near  the  epicenters. 


£.  Attenuation  at  the  Alzerlan  Test  Site 


There  have  been  several  studies  in  which  the  amplitudes  of  the  seismic  waves 
emitted  from  explosions  at  the  French  test  site  in  Algeria  were  measured  to  esti- 
mate the  seismic  yields  from  these  explosions.  In  these  studies  the  tectonics  of 
the  Algerian  site  were  assumed  to  be  typical  of  a continental  shield.  In  an  ex- 
tensive review  on  Intraplate  seismicity,  Sykes  (1978)  points  out  that  the  region 
of  the  French  test  site  has,  contrary  to  previous  belief,  not  been  tectonically 
stable  since  the  Precambrlan.  The  region  was  reactivated  during  the  separation  of 
Africa  from  the  Americas,  and  has  experienced  volcanlsm  during  the  past  25  million 
years.  Such  regions  of  reactivation  tend  to  follow  old  zones  of  weakness  along 
the  edges  of  cratons.  The  French  test  site  is,  therefore,  located  in  a zone  that 
has  experienced  tectonic  activity,  including  volcanlsm  and  uplift,  during  the  past  25 
million  years.  Seismic  waves  leaving  this  region  may  be  severely  attenuated  in  ampli- 


tude. This  situation  is  quite  similar  to  the  Nevada  Test  Site,  a region  that  has  also 


i.. 


experienced  late  Tertiary  and  Cenozoic  volcanism.  Hence,  amplitude-yield  curves 
for  the  Algerian  test-site  probably  cannot  be  directly  applied  to  estimating 
yield  for  sites  that  have  been  tectonically  stable  for  hundreds  of  millions  of 


years. 


II.  RELATED  STUDIES 


The  following  section  summarizes  research  that,  although  not  supported 
directly  by  the  contract.  Is  of  Interest  to  the  subject  of  the  contract. 
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A.  Tectonics  of  Central  Asia 


The  tectonics  of  central  Asia  are  complicated.  In  contrast  to  the  relatively 
simple  structures  observed  at  subductlon  zones,  those  structures  that  were  pro- 
duced In  response  to  the  collision  between  continental  India  and  Eurasia  are  not 
as  well  understood.  In  this  section  we  describe  several  projects  carried  out  at 
L-DGO  that  concern  central  Aslan  tectonics. 

1.  Seismic  activity  recorded  by  a network  of  telemetered  stations  In  northern 
Pakistan  has  recently  revealed  a detailed  picture  on  the  structures  and  processes 
currently  affect  ng  the  Hazara  region.  By  unraveling  the  complicated  situation 
that  exists  there.  It  may  now  be  possible  to  understand  the  tectonics  of  other 
regions  in  central  Asia  that  are  characterized  by  a similar  environment.  Such 
regions  may  Include  the  eastern  termination  of  the  Himalayan  structures  where  sev- 
eral anomalous  earthquakes  have  occurred. 

The  Hazara  region  of  northern  Pakistan  Is  an  arcuate  belt  of  foreland  folds 
and  thrusts  that  extends  from  the  western  Himalayan  syntaxls  to  the  Sulalman 
Range  (Figure  6A,  B) . A model  based  on  the  mapped  distribution  of  mlcroselsmlclty, 
composite  fault-plane  solutions,  and  geological  evidence  was  constructed  to  explain 
the  recent  tectonic  processes  of  this  region  (Seeber  and  Jacob,  1977;  Armbruster 
et  al. . 1978;  Seeber  and  Armbruster,  1979).  The  Indus-Kohlstan  Seismic  Zone  and 
the  Hazara  Lower  Seismic  Zone  are  major  zones  of  displacement  within  the  basement 
(Figure  7);  the  former  Is  the  more  active  of  the  two.  These  zones  appear  to  be 
extensions  of  the  Himalayas  beyond  the  western  Himalayan  syntaxls.  The  Hazara 
Lower  Seismic  Zone  seems  to  merge  upward  and  to  the  south  with  a shallow-dipping 
surface  of  decollement  that  decouples  the  basement  from  the  overlying  sediments 
and  metasediments  (which  have  a thickness  of  between  10  to  30  km)  (Figure  8) . 

The  decollement  surface  Is  probably  lubricated  by  Infracambrlan  salt. 
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2.  On  a braoder  scale,  the  historical  and  modern  seismicity  (Figures  9 and  10  ) 
for  a portion  of  south  central  Asia  has  been  compiled  and  critically  reviewed  to 
reveal  the  large-scale  selsmotectonic  trends  (Qulttmeyer  and  Jacob,  1979;  Qulttmeyer 
et  al. , 1979).  It  is  difficult  to  associate  individual  earthquakes  with  specific 
mapped  faults  for  the  following  reasons:  (1)  areas  are  mapped  only  in  a recon- 
naissance fashion,  (2)  the  locations  of  the  earthquakes  are  of  limited  precision, 
and  (3)  the  density  of  faults  is  sometimes  high.  In  some  cases,  however,  seismi- 
city can  be  associated  with  faults  or  zones  of  faulting.  This  seismicity  includes 
large  earthquakes  and  their  associated  aftershocks,  and  well  defined  llneatlons  of 
earthquakes  with  small  to  moderate  magnitudes. 

3.  The  Makran  region  of  southern  Pakistan  has  been  the  subject  of  two  moderately 
detailed  studies.  Jacob  and  Qulttmeyer  (1979)  showed  that  volcanic,  geologic  and 
seismic  evidences  are  consistent  with  the  interpretation  of  the  Makran  region  as 
a zone  of  active  subductlon.  If  rigid  plate  assumptions  are  made  and  the  poles  and 
rotation  rates  of  McKenzie  (1972),  McKenzie  and  Slater  (1971)  and  Minster  et  al. 
(1974)  are  used,  it  is  found  that  oceanic  portions  of  the  Arabian  plate  subduct 
northwards  towards  Eurasia  at  a rate  of  about  5 cm/yr.  Most  of  the  tectonic  fea- 
tures associated  with  the  arc-trench  system  are  extreme.  For  instance: 

(1)  the  arc-trench  gap  measures  500  + 100  km  in  width,  (2)  the  Benloff  zone  is  only 
weakly  developed  to  a depth  of  80  km,  (3)  the  volcanic  chain  shows  wide  spacing 
(>  100  km)  between  major  centers  of  activity,  and  (4)  a large  portion  of  the  accre- 
tionary prism  is  subaerlally  exposed  (see  Figures  11  and  12  ) . Two  focal  mechanisms 
in  the  deepest  part  of  the  documented  Benloff  zone  exhibit  down-dip  tension.  One 
great  earthquake  (M^  * 8)  in  1945  ruptured  a section  of  the  plate  boundary. 

4.  Seismic  activity  prior  to  the  1945  earthquake  along  the  Makran  coast  of 
southern  Pakistan  was  consistent  with  patterns  identified  by  Kelleher  and  Savino 
(1975)  and  Mogi  (1969)  before  other  large  earthquakes  (Qulttmeyer,  1979).  Although 


the  data  are  sparce,  the  seismic  activity  inland  (probably  within  the  subducted 
slab)  for  earthquakes  with  magnitude  6 to  7 appears  to  be  higher  for  the  period 
20  to  10  years  before  the  great  coastal  earthquake  than  in  any  10  year  period 
after  it.  In  addition,  activity  occurred  along  the  coast  only  in  the  vicinity 
of  the  impending  1945  epicenter  (Figure  13A-D) . The  region  that  was  selsmlcally 
active,  both  before  and  after  the  1945  shock,  is  now  relatively  quiet. 


B.  Tectonics  of  the  Aleutian  Seismic  Zone 


A number  of  studies  concerning  the  Aleutian  seismic  zone  are  currently  being 
pursued  at  L-DGO.  These  studies  Increase  our  knowledge  of  the  configuration  of 
a typical  subductlon  zone  and  some  of  the  Interesting  earthquakes  that  have  occur- 
red in  It.  The  results  of  these  studies  may  be  applicable  to  understanding  other 
subductlon  zones  such  as  the  Kuriles  and  Kamchatka. 

1.  A seismic  network  has  been  operated  In  the  Shumagln  Islands  and  adjacent 
Alaska  Peninsula,  Alaska  since  July  1973.  Hypocentral  cross-sections  based  on 
the  first  four  years'  data  from  this  network  show  a well-defined  Benloff  zone  that 
Is  about  10  km  thick  and  extends  to  about  180  km  depth  (Davies  and  House,  1978) . 

Focal  mechanism  and  strong-motion  accelerograph  data  for  an  m^  * 6.0  earthquake 
at  40  km  depth  Indicate  a slip  plane  subparallel  to  the  Benloff  zone.  This  event 
Is  Interpreted  as  being  directly  caused  by  the  underthrusting  of  the  Pacific 
lithosphere  beneath  the  Shumagln  Islands.  Comparison  of  the  Shumagln  hypocentral 
cross-section  with  cross-sections  from  other  regions  of  the  Aleutian  arc  show  that 
at  depths  greater  than  40  km  (1)  the  Benloff  zones  are  nearly  congruent  and  (2) 
their  positions  relative  to  the  volcano  line  are  nearly  the  same  along  the  entire 
length  of  the  arc.  This  comparison  also  reveals  that  the  eastward  widening  of  the 
volcano-trench  separation  occurs  predominantly  In  the  shallowly  dipping  extension 
of  tha  Benloff  zone  from  Its  position  at  40  km  depth  to  where  It  outcrops  near  the 
trench;  this  extension  Is  herein  termed  the  isaln  thrust  zone.  Rupture  zones  of  the 
1964,  1965,  and  possibly,  the  1938  great  earthquakes  along  the  Aleutian  arc  are 
confined  to  the  main  thrust  zone.  This  observation  and  the  very  low  level  of  seismic 
activity  at  the  present  time  In  the  main  thrust  zone  compared  to  that  In  the  Benloff 
zone  below  40  km  depth  suggest  that  seismic  activity  In  the  main  thrust  zone  Is 
dominated  by  great  earthquakes  and  their  aftershocks,  followed  by  long  Intervals 
of  quiescence.  The  Benloff  zone  below  40  km  depth  shows  more  continual  seismic 
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activity.  For  the  Aleutian  arc  there  exists  a weak  correlation  between  the  width  of 


the  main  thrust  zone  and  the  rupture  length  and  magnitude  (Mw)  of  great  earthquakes. 
2.  Earthquakes  with  high  stress-drop  generate  an  unusually  large  amount  of  high- 
frequency  seismic  waves.  They  also  have  a large  m^  compared  to  M^.  Since  these 
seismic  characteristics  are  similar  to  those  of  nuclear  explosions,  close  examina- 
tion of  such  high  stress-drop  earthquakes  may  further  our  understanding  oii  the 
limitation  of  this  discrimination  method. 

Two  moderate-size  earthquakes  (mj^  ■ 5.8,  6.0),  with  high  estimated  stress  drops, 
occurred  on  April  6,  1974,  beneath  the  Shumagin  Islands  (Alaska)  within  a local  net- 
work of  short-period  seismograph  stations.  These  earthquakes  triggered  a strong- 
motion  accelerograph  (SMA)  at  Sand  Point  (SDF)  (Figure  14  ) . 

P arrivals  for  the  two  mainshocks  and  29  aftershocks  were  picked  with  a pre- 
cision of  0.05  sec  on  seismograms  recorded  by  5 local  stations.  The  locations  of 
these  earthquakes  are  presented  in  a cross-sectional  view  in  Figure  15.  The  hypo- 
centers  define  a plsme  dipping  at  about  30°  to  the  NW  (Figure  15E) . Figure ri5A  is  a 
view  normal  to  the  local  Wadati-Benioff  zone.  This  figure  shows  Chat  the  after- 
shock zones  for  the  two  mainshocks,  although  close  together,  are  relatively  dis- 
tinct. Linear  dimension  (diameter)  of  the  aftershock  zone  is  about  3-4  km  for  the 


m^  ■ 5.8  event 

Long  period  data  from  the  WWSSN  were  used  to  determine  the  focal  mechanism  of 


and  about  4-5  km  for  the 


■ 6.0  event. 


Che  0356,  m^  ■ 6.0  earthquake  (Figurel6). 
ficient  to  constrain  either  nodal  plane. 


P-wave  first  motions  alone  are  Insuf- 
The  addition  of  S-wave  polarizations  re- 


fines the  mechanism  determination  considerably  and  places  limits  on  its  uncertainty. 
The  tectonic  setting  and  the  trend  of  the  aftershocks  favor  the  choice  of  the  NW 


dipping  nodal  plane  as  the  fault  plane.  Focal  parameters  for  the. fault  plane  are: 
strike  ■ 254*  + 15°,  dip  ■ 30°  + 5°,  and  rake  of  slip  vector  ■ 90*  + 15*. 

Analysis  of  the  SMA  records  was  carried  out  using  a circular-fault  version  of  a 
qw.si-dynamic  model  (Boatwright,  1978b,  1979).  A description  of  this  technique. 


summarized  from  Boatwright  (1978b,  1978),  follows: 

The  S-wave  signals  were  digitized,  corrected  for  Instrument-response,  and  ro- 
tated to  and  components.  Since  the  amplitudes  were  less  than  lOZ  of 
those  for  S^,  the  analysis  was  restricted  to  the  component  only. 

Accelerograms,  corrected  for  the  effect  of  free  surface,  are  shown  In  Figures  17 

2 

and  18,  as  well  as  respective  velocity  and  V traces.  Note  the  remarkable  similarity 
2 

of  the  V plots  In  pulse  shapes  and  amplitudes  for  the  two  events.  The  events 

probably  share  a similar  rupture  geometry  and  velocity. 

Detailed  source  models,  derived  from  the  circular  version  of  the  quasl-dynamlc 

model,  have  an  "elliptical"  or  self-similar  slip  distribution  and  causal  healing. 

-2 

The  displacement  spectra  for  these  two  events  fall  off  more  steeply  than  u . The 
sharp  falloff  Indicates  a gradual  stopping  of  the  rupture.  This  feature  was  also 
Incorporated  Into  the  model.  Final  model  parameters  are  listed  In  Table  2. 

The  high  stress-drops  for  these  two  earthquakes  may  Indicate  either  a regional  stress 
bulld-up  prior  to  a large  earthquake  or  the  release  of  a localized  high  strength 
region  (asperity)  that  was  unbroken  by  the  most  recent  (1938)  major  earthquake. 

The  calculated  stresses  are  among  the  highest  computed  for  any  seismic  gap  In  the 
world.  Archambeau  (preprint)  has  also  obtained  a high  stress-drop  estimate  for 
several  recent  earthquakes  In  the  Shumagln  region. 

3.  Seismicity  located  using  the  Adek  seismograph  network  shows  the  Benloff 
zone  below  a depth  of  100  km  consists  of  two  thin  zones  of  earthquakes  about  25  km 
apart  that  merge  at  a depth  of  175  km.  Focal  mechanisms  In  the  upper  zone  are 
consistently  down-dip  compression,  while  those  of  the  lower  zone  are  down-dip  ten- 
sion. An  elastic  plastic  model  of  the  lithosphere  predicts  that  these  two  zones 
reflect  the  stresses  In  the  elastic  core  of  the  lithosphere  due  to  unbending 
(Engdahl  and  Scholz,  1977). 


C.  North-Atlantlc  Intraplate  Earthquake  of  24  March  1978;  Preliminary  Results 


On  March  24,  1978,  an  Intraplate  earthquake  of  magnitude  6 occurred  about 
350  km  southwest  of  Bermuda,  In  seafloor  of  Mesozoic  age  (near  magnetic  anomaly 
M4:  'V'  117  ay  B.F.)>  Teleselsmlc  records  of  this  shock  were  analyzed  to  obtain  a 
fault-plane  solution.  An  aftershock  survey  was  conducted  from  June  17-24,  1978, 
to  study  the  settings  and  possible  causes  of  this  Intraplate  earthquake.  Since 
several  test  sites  also  have  Intraplate  environments  (e.g.  - MTS,  Algeria,  India 
and  several  In  the  U.S.S.R.),  It  Is  Important  to  understand  the  causes  and  charac- 
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terlstlcs  of  Intraplate  earthquakes  so  that  we  may  discriminate  them  from  nuclear 
tests  in  Intraplate  environments. 

The  aftershock  survey,  conducted  In  cooperation  with  Woods  Hole  Oceanographic 
Institution,  deployed  sonobuoy  and  ocean  bottom  hydrophone  (OBH)  arrays.  Con- 
flicting locations  reported  for  the  malnshock  resulted  In  the  misalignment  of  the 
array  away  from  the  center  of  aftershock  activity.  Nevertheless,  numerous  events 
were  recorded.  Preliminary  results  place  the  cluster  of  aftershock  activity  at 
29*50'N,  67*20'W. 

The  main  shock  was  of  thrust  type.  Data  at  hand  are,  however,  lnsufflcie:.t  to 
constrain  the  orientation  of  the  nodal  planes.  The  Investigation  of  this  earth- 
quake Is  continuing. 

This  series  of  shocks  has  an  excellent  set  of  high  frequency  observations  at 
Intermediate  distances  (about  15*  to  25*)  which  are  of  great  Importance  to  the 
understanding  of  the  detection-discrimination  levels  that  a network  In  the  U.S.S.R. 
could  provide. 


Two  nearby  earthquakes  seaward  of  the  Tonga  trench,  one  normal-faulting  and 
one  thrust-faulting,  were  studied  In  detail  using  comparative",  surface  wave  tech- 
niques. The  techniques  require  the  epicenters  of  two  earthquakes  to  be  close  to 
each  other  so  that  the  effects  of  propagation  and  station  response  for  common 
stations  would  be  the  same  for  both  events.  By  looking  only  at  relative  ampli- 
tudes of  the  two  events  at  each  station,  much  of  the  scatter  can  be  removed.  In 
addition,  since  the  propagation  paths  are  similar,  the  relative  phases  of  the  Ray- 
leigh waves  at  each  station  can  be  compared.  The  observed  amplitudes,  corrected 
for  the  effects  of  propagation  and  station  response,  and  observed  focal  phases  were 
fitted  with  theoretical  patterns  to  obtain  the  best-fitting  focal  mechanism  solu- 
tions. The  focal  depth  for  the  thrust-faulting  event  was  found  to  be  49  km  below 
sea-level,  consistent  with  the  focal  depth  estimated  from  pP-F  delay  time.  The 
normal-faulting  event  had  a depth  of  14  km. 

These  are  the  first  precise  depth  determinations  for  earthquakes  which  are 
associated  with  the  bending  of  oceanic  lithosphere  before  subductlon.  The  pattern 
of  horizontal  devlatorlc  tension  near  the  surface  and  compression  within  the  In- 
terior of  the  lithosphere  Is  consistent  with  the  stresses  predicted  by  an  elastic 
plate  model,  although  the  depth  of  the  thrust  event  suggests  chat  either  the 


elastic  layer  begins  near  the  base  of  Che  crust  or  chat  the  rheological  model  should 
be  modified  to  Include  a plastic  layer  of  finite  strength  (Chen  and  Forsyth,  1978) . 


Important  aspect  of  a nuclear  discrimination  program. 
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TABLE  2 

Derived  Source  Parameters 


0153-event 
Radius  - a ~ 1.2  km 

24 

Moment  -M  “3.5+.8xl0  dyne • cm 
o — 

Stress  drop  - Ao  ■ 890  bars  Range  600-1100  bars 

Effective  Stress  - t ■ 1040  + 350  bars 

® 20 

Radiated  Energy  -E  *8. 7+3. Ox  10  dyne* cm 

Apparent  Stress  “ i bars 


0356-event 

Radius  - a ■ 1-65  km 

24 

Moment  -M  ■ 6.7  +1.5  x 10  dyne* cm 
o — 

Stress  Drop  - Ao  ~ 650  bars  Range  350-800  bars 

Effective  Stress  - t ■ 780  + 250  bars 

® 20 

Radiated  Energy  -E  ■12. 4+4. Ox  10  dyne* cm 
Apparent  Stress  ” ■ 120  + 50  bars 
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FIGURE  CAPTIONS  i 

Figure  1.  a.  Percentage  of  energy  transmitted  as  a function  of  angle  in  poly- 
styrene. 

b.  Percentage  of  energy  reflected  as  a function  of  incident  angle. 

Figure  2.  Group  velocities  for  Rayleigh  and  Love  waves  across  the  Tibetan  Plateau. 

a.  The  earthquakes  (filled  circles)  and  WWSSN  stations  (open  circles) 
used  in  the  study.  The  Tibetan  Plateau  is  strlppled;  the  Indo-Gungetic 

basin  is  hatched.  , 

4 

b.  Group  velocities  for  the  paths  in  (a).  The  smallest  dot  represents 
one  data  point;  the  largest  dot  represents  three  or  more  data  points 
superimposed. 

c.  Averaged  group  velocities  for  Rayleigh  and  Love  waves.  The  bars 
indicate  twice  the  RMS  value  of  the  residuals.  Previously  published 
Rayleigh-wave  data  by  Gupta  and  Naraln  (1967),  Santo  and  Sato  (1966) 
and  Chen  and  Molnar  (1975)  are  also  shown. 

Figure  3.  a.  Acceptable  shear  wave  models  (TP-2,  TP-3,  and  TP-4)  for  the  Tibetan  ^ 

Plateau.  Note  the  low-velocity  zone  at  depths  of  about  25  to  40  km.  j 

Also  shown  for  comparison  is  the  Canadian  Shield  model  CANSD  (Brune  and  # 

Dorman,  1963) . 

b.  Vertical  cross  section  along  85*’E  as  shown  in  Figure  3a  by  the  solid 
line  X-X' . TF-4  is  from  this  study.  Other  seismic  models  compiled  are: 

I 

(1)  Bhattacharya  (1971),  (2)  Gangetlc  model  of  Chun  and  Yoshii  (1977),  | 

(3)  Chen  and  Molnar  (1975)  and  (4)  Gupta  and  Narain  (1967).  Depth  to 
the  "Moho"  is  from  a Chinese  tectonic  map  (Tectonic  Map  Compiling  Group,  f 

1974  - solid  line),  and  deduced  from  gravity  data  (Kono,  1974  dashed 


line) . 


Figure  4 
Figure  5 


Figure  6. 


Figure  7. 
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. See  Figure.  , ! 

. a.  Approximate  location  of  earthquakes  studied  (filled  circle)  and 
the  WWSSN  stations  used  in  surface  wave  study  of  shocks  from  mid- 
Atlantic  Ridge. 

b.  - e.  The  amplitude  spectral  densities  plotted  as  a function  of  azi- 
muth for  events  1 to  4,  respectively.  The  spectral  densities  are  given 
at  the  following  periods:  event  1,  35.7  sec,  event  2,  35.7  sec,  event  I 

3,  35.7  sec,  event  4,  65.0  sec.  Event  1 (Figure  5B)  shows  a two-lobed  j 

radiation  pattern  typical  of  shallow  normal-faulting  earthquakes. 

Event  4E  shows  a four-lobed  radiation  pattern,  signifying  strike-slip 

faulting.  Radiation  patterns  for  events  2 and  3 are  ambiguous,  and  | 

cannot  resolve  the  fault  type. 

a.  Regional  setting  of  the  Hazara  region  in  northern  Pakistan.  Areas 
other  than  Quaternary  basins  are  dotted. 

b.  Seismic  stations  (triangles)  comprising  the  Tarbela  Dam  (TD)  and  | 

j 

Chashma  Power  Plant  (near  CR)  networks.  Structural  features  are  shown  | 

for  reference  and  are  the  same  as  shown  in  Figure  7.  | . 


Six  months  of  eplcentral  data  from  the  Tarbela  and  Chashma  networks. 

Only  earthquakes  of  magnitude  0.5  or  greater  are  shown.  The  seismicity 
plotted  is  representative  of  the  longer  term  activity  as  recorded  by 

the  networks.  The  concentration  of  activity  in  the  upper  right  that  has 
a northwest-southeast  trend  is  the  Indus  Kohlstan  Seismic  Zone  (IKSZ). 

The  Hazara  Lower  Seismic  Zone  parallels  the  IKSZ  about  80  km  to  the  south- 
west. Geologic  features  are  shown  in  thinner  lines;  faults  identified 
from  seismic  data  are  shown  in  thicker  lines.  Line  X-X'  represents  the 
projection  line  of  the  cross-section  shown  in  Figure  8. 


f 

I 

I • 

f 

I 

I-. 

J ' 

?■'■■■ 


i 


Figure  8.  A vertical  projection  of  the  seismicity  along  a 20  km  wide  northeast- 
trending strip  across  the  Tarbela  network.  All  good  hypocenter  locations 
between  August  1973  and  August  1976  are  included.  Solid  triangles  are 
seismic  stations.  The  average  topographic  relief  Is  Indicated.  Composite 
fault-plane  solutions  are  projected  on  hemispheres  convex  to  the  south- 
east (towards  the  viewer) . Solution  12  Is  for  earthquakes  west  of  the 
HLSZ  west  of  the  Indus  River  (not  shown  In  the  section) . Faults  are 
drawn  as  suggested  by  h3rpocenter  alignments;  directions  of  motion  are 
from  the  composite  solutions.  The  thrusting  motion  Indicated  for  the 
low-angle  fault  Is  conjectural. 

Figure  9.  Map  of  maximum  documented  Intensity  (Modified  Mercalll  Scale)  at  any 
given  location.  Data  for  the  time  period  25  A.D.  to  1972  are  shown; 
however,  mapped  portions  of  the  U.S.S.R.  and  China  are  not  considered, 
and  data  from  Iran  are  Insufficient  for  quantitative  description. 
Isoselsmal  lines  (dotted  where  Inferred)  are  plotted  for  some  of  the 
larger  events.  The  year  of  occurrence  for  each  such  large  event  Is 
Indicated.  The  Intensity  value  associated  with  a given  Isoselsmal  line 
Is  Indicated  In  the  box  near  each  date.  The  first  value  given  Is  for 
the  Innermost  isoselsmal  line,  etc.  A few  locations  for  which  a docu- 
mented Intensity  Is  known  are  not  plotted  so  that  Isoselsmal  lines  will 
be  more  clearly  visible.  The  open  triangles  represent  some  major  cities. 


Figure  10.  Eplcentral  map  of  modern  seismicity  for  the  region  studied.  The  filled 
squares  are  proportional  to  magnitude  and  represent  earthquakes  recorded 


at  teleselsmlc  distances  from  January,  1914  through  April,  1975.  The 
symbol  labelled  UD  In  the  magnitude  scale  denotes  events  with  undeter- 
mined magnitude.  Events  from  1914  to  1964  are  relocated.  No  distinction 
Is  made  between  surface-wave  and  body-wave  magnitudes  In  this  figure. 

Four  large  earthquakes  (M^  >7)  that  occurred  from  1905  to  1914  are 
represented  by  open  circles. 

Figure  11.  Tectonic  sketch  map  of  the  Makran  region  and  surrounding  vicinities. 

The  abyssal  part  of  the  oceanic  portion  of  the  Arabian  plate  (In  the 
Gulf  of  Oman),  that  Is,  the  area  with  water  depths  larger  than  3000  m. 

Is  shown  In  dark  diagonal  hatching.  Note  the  distribution  of  ophlollte 
belts  and  "colored  melanges"  on  land  (solid  dark  areas).  In  some  cases, 
they  may  Indicate  the  approximate  locations  of  probable  former  continent/ 
ocean  boundaries  In  the  late  Mesozoic /early  Cenozolc.  Also  note  that  the 
pattern  of  present  drainage  divides  (light  ring  patterns)  mimics  approxi- 
mately the  plate  boundaries  and  the  triple  junction,  probably  as  they 
may  have  existed  some  30  (or  more)  m.y.  ago.  Deep  sediment-filled 
Interior  basins,  probably  with  recent  histories  of  subsidence,  are  shown 
In  light  shading.  Of  these,  the  Jaz-Murlan  and  the  Hamun-l-Mashkel 
depression  are  Interpreted  by  Farhoudl  and  Karlg  (1977)  as  forearc 
basins;  they  lie  seaward  of  the  volcanic  arc.  Note  the  linear  arrange- 
ment of  the  volcanic  edifices  In  the  approximately  450  km  long  and  less 
than  100  km  wide  volcanic  arc.  The  three  major  Quaternary  volcanic 
centers  - the  Koh-l-Sultan  (elevation  2,510  m),  Taftan  (4,100  m),  and 
Bazman  (3,490  m)  - are  spaced  more  than  100  km  apart,  but  connected  by 
numerous  small  volcanic  plugs  and  edifices.  The  volcanic  arc  strikes 
approximately  ENE  and  cuts  diagonally  across  numerous  tectonic  surface 
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features.  The  carbonatlte  volcano  Koh-l-Khan  Nasln  In  the  Helmand 
Basin  of  Afghanistan  (upper  center)  Is  probably  not  part  of  the 
Baluchistan  volcanic  arc  which  produces  mostly  rocks  of  andesitic, 
calc-alkallne  composition. 

Figure  12.  Cross  section  XX'  (along  bO^E  from  24”  to  SCN)  through  the  western 

Makran  region  showing  the  earthquake  hypocenters  and  Inferred  dipping 
Benlof f zone  (bottom) , and  topography  and  some  surface  tectonic  fea- 
tures (top)  along  the  same  profile.  Bottom;  Circles  represent  events 
up  to  200  km  to  the  east  of  the  section  line,  triangles  up  to  200  km 
to  the  west.  Filled  symbols  represent  events  for  which  the  depth  Is 
constrained  by  at  least  one  depth  phase;  open  symbols  represent  events 
for  which  the  depth  Is  determined  by  minimizing  the  residuals  of  first 
P arrivals  only.  Average  depths  determined  Independently  by  these  two 
methods  usually  differ  by  no  more  than  10  km  In  cases  when  both  are  avail- 
able In  this  region.  The  events  labelled  A and  B are  those  for  which 
fault  plane  solutions  were  determined.  The  arrows  at  these  two  hypo- 
centers  show  the  plunges  of  the  T axes  for  these  two  events,  taking  Into 
account  the  2X  vertical  exaggeration.  The  shaded  band  Is  the  Inferred 
upper  boundary  of  the  descending  oceanic  lithosphere  belonging  to  the 
Arabian  plate.  Top;  Major  tectonic  features  and  subdivisions  of  the 
trench-volcano  gap  In  the  western  Makran.  Note  the  coincidence  In  the 
portion  of  trench-slope  break  at  the  southern  boundary  of  the  subsiding 
forearc  basin  (Jaz  Murlan  depression)  with  the  sudden  deepening  of  the 
Benlof f zone  from  crustal  (<  40  km)  to  sub-crustal  depths.  The  develop- 
ment of  these  forearc  basins  n»y  be  closely  related  to  a sudden  steepening 
In  the  dip  of  the  descending  oceanic  lithosphere  beneath  It  (Jacob  ^ , 

1977).  Note  the  different  vertical  exaggerations  in  the  top  and  bottom 
of  the  Figure. 


Figure  13.  Seismicity  of  the  Makran  region  during  different  time  periods.  The 
magnitude-quality  scale  (note:  there  Is  no  differentiation  among  A 
B,  and  C quality  events)  Is  shown  In  part  A;  the  distance  scale  Is 
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shown  In  part  B.  Events  for  which  a magnitude  could  not  be  determined 
are  plotted  as  If  they  had  a magnitude  of  < 4.9.  The  epicenter  of  the 
great  earthquake  In  1945  Is  represented  by  a cross  enclosed  In  a circle 
In  parts  A,  C,  and  D.  The  major  centers  of  Quaternary  volcanlsm  are 
Indicated  by  filled  circles  surrounded  by  radiating  line  segments. 

A.  Seismic  activity  before  the  1945  earthquake.  Activity  Is  located 
north  and  northwest  of  the  1945  rupture  zone  and  in  the  vicinity  of 
the  1945  epicenter.  The  region  west  of  this  activity  is  relatively 
quiet.  The  two  earthquakes  In  the  southeast  corner  of  the  map  are 
associated  with  the  Murray  ridge. 

B.  The  great  1945  earthquake  and  its  long-term  aftershock  activity.  The 
largest  symbol  (filled  circle)  represents  the  epicenter  of  the  1945  event. 
The  long-term  aftershock  activity  suggests  a rupture  length  of  100  to  200 
km  to  the  east  of  the  main  shock  epicenter.  The  region  northwest  of  the 
rupture  zone  Is  now  relatively  quiet;  the  region  west  of  the  rupture  zone 
remains  quiet. 

C.  Seismicity  after  the  decline  In  long-term  aftershock  activity.  During 
this  period  of  time  the  entire  Makran  region  exhibits  a low  level  of 
activity. 

D.  Recent  seismicity  of  the  Makran  region.  An  Increased  number  of  earth- 
quakes are  detected  because  of  the  establishment  of  the  WWSSM.  Most  of 
these  events  would  not  have  been  detected  during  an'  earlier  time  period. 
Activity  along  the  coast,  west  of  the  1945  epicenter  produces  a "donut" 
pattern.  One  earthquake  of  moderate  magnitude,  located  to  the  northwest 
of  the  donut  pattern,  may  be  analogous  to  the  activity  that  preceded  the 
1945  event. 
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Figure  4a*  Normalized  travd-time  re- 
siduals at  individual  stations  in  northeast- 
etn  United  States  and  adjacent  parts  of 
Canada  for  Cannikin,  Milrow,  and 
Longshot  explosions  in  Aleutians.  Indi- 
vidual readings  and  standard  deviations  are 
listed  in  Table  1 (see  footnote  1).  Dashed 
line  boundary  between  rocks  of  Frc- 
cambrian  age  to  north  and  younger  rocks  to 
south.  Note  zone  centered  in  northetn  New 
York  and  southern  Quebec  where  all  re- 
siduals are  unusually  negative  by  -0.6  to 
— 1.5  s.  Larger  numbers  indicate  that  re- 
sidual was  determined  from  two  or  more 
explosions  with  standard  deviation  less 
than  0.35  s.  Smaller  symbols  denote 
readings  from  either  single  explosions  or 
averages  with  larger  standard  deviations. 


Figure  4b.  Normalized  travel-time  re- 
siduals for  explosions  at  Novaya  Zemlya 
USSR  (in  parentheses)  and  Nevada  Test 
Site.  Size  of  symbols  same  as  Figure  4e 
Data  in  Table  1 (see  footnote  1). 
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Figure  lA.  Location  nap  of  the  Shumagin  Islands  and  April  6,  1974  earthquakes  (solid  dot  in  circle).  Shumagin 

network  seisnic  stations  are  indicated  by  triangles  and  three  letter  station  circles.  For  comparison 
the  (mislocated)  PDE  epicenter  is  also  shown  (solid  square  in  circle). 
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APRIL  6.  19  74  0356 

54.91  N I60.28W  40  KM 

Figure  16.  Focal  mechanisa  of  the  0356  hra  event  on  April  S,  1974.  Flrat  motlona  and  S-wave 
polarlzatlona  are  projected  onto  a lover  hemlaphere,  equal  area  net.  WWSSN  long 
period  flrat  motlona  are  clrclea,  filled  are  compreaalona,  empty  dllatatlona; 
aaaller  ayabola  ere  leaa  reliable  data.  Local  netvork  abort  period  flrat  motlona 
are  aquarea;  filled  are  compreaalona;  empty,  dllatatlona.  S-vava  polarlzatlona  are 
llnea  on  figure,  leaa  reliable  determlnatlona  are  only  half  the  length. 
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Figure  17.  S componenC  of  S wave  for  0153  event,  April  6,  1974.  Bottom:  Instrument 
and  free  surface  corrected  acceleration.  Center:  Velocity  trace  obtained 
from  Integrating  baseline  corrected  acceleration.  Top:  Velocity-squared 
trace;  dashed  line  is  velocity-squared  trace  of  final  model. 
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ts.  supplementary  notes 


IS.  KEY  WORDS  fConlInu#  on  fOYoeoo  otdo  II  nocot oory  end  Idontliy  by  block  number) 

Reflection/ transmission  coefficients,  numerical  methods,  group  velocities,  Ray- 
leigh waves.  Love  waves,  surface  waves,  crustal  velocities,  Tibet,  P-wave  resi- 
duals, suture  zones.  New  York  State,  Vermont,  Quebec,  focal  mechanism,  Mid- 
Atlantic  ridge,  Algeria,  magnitude-yield  relation,  intraplate  seismicity.  Cen- 
tral Asia,  Pakistan,  composite  fault-plane  solutions,  mlcroselsmlcity , Himalayas 

20.  ABSTRACT  (Conilnuo  on  rovoroo  tide  II  nocoooory  and  Idontliy  by  block  numbor) 

The  results  of  a number  of  projects  are  summarized: 

An  approximate  method  for  calculating  the  reflection  and  transmission 
coefficients  of  obliquely  incident  Rayleigh  waves  across  two  welded 

quarter-spaces  Is  outlined. 

/ 

Group  velocities  of  fundamental  mode  Rayleigh  and  Love  waves,  with 
paths  traversing  the  Tibetan  Plateau,  were  used  to  determine  that  the 
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Block  19  (con't.) : 

Hazara  arc,  historical  seismicity,  Makran,  subductlon,  volcanlsm,  Benloff  zone, 
seismic  gaps,  Aleutians,  stress-drop,  strong-motion  accelerograms,  lithosphere, 
aftershocks. 


Block  20  (con't.): 

^Tibetan  Plateau  has  a 70  km  thick  crust,  with  a low  velocity  zone  at  an 
Intermediate  level  within  the  crust. 

y 

A zone  of  large,  negative  residuals  for  P-wave  arrivals  exists  In  north- 
ern  New  York  State. 

^ — 4*  Surface-wave  amplitude  spectra  constrained  by  body-wave  first  motions 
arc  sufficient  to  determine  focal  mechanisms  for  some  moderate-sized  earth- 
rtuakes  along  the  Mid-Atlantic  ridge,  'n 


S.  The  French  test  site  In  Algeria  has,  contrary  to  previous  belief,  not 
been  tectonically  stable  since  the  Precambrlan.  Thus,  seismic  waves  leaving 
this  region  may  be  more  attenuated  than  was  previously  suspected.  Yield 
estimations  for  the  French  test  site  will  have  to  be  re-evaluated. 

^ fty’^lcro seismicity  and  geological  evidences  were  used  to  model  the  current 


tectonic  processes  of  the  Hazara  district  In  northern  Pakistan. 


^ wa 


■fv  The  historical  and  modern  seismicity  of  a portion  of  south-central  Asia 
was  compiled  and  critically  reviewed  to  delineate  large-scale  selsmotectonlc 
trends.  a^cSL 


Seismic,  volcanic,  and  plate  tectonic  data  Indicate  the  Makran  region  of 
' southern  Pakistan  Is  an  active  subductlon  zone.  ^ 

9.  Seismicity  and  tectonics  of  the  Aleutians,  Alaska  were  studied.  Large 
stress-drops  (“v  500  bars)  for  two  earthquakes  were  determined  from  strong- 
motion  accelerograph  records. 

10.  A focal  mechanism  study  and  aftershocks  survey  were  undertaken  for  a 
large  Intraplate  earthquake  350  km  southwest  of  Bermt^da. 

11.  A comparative  surface-wave  technique  was  used  to  precisely  determine  the 
depths  of  two  earthquakes  that  occurred  seaward  of  the  Tonga  trench.  Im- 
plied stresses  for  these  earthquakes  are  consistent  with  bending  of  an 
elastic  plate  model. 
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